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Abstract 

An analytical study was performed as part of the NASA 
Lewis support of a National Transportation Safety Board 
(NTSB) aircraft accident investigation. The study was 
focused on the performance degradation associated with ice 
contamination on the wing of a commercial turbo-prop-pow- 
ered aircraft. Based upon the results of an earlier numerical 

study conducted by the authors 111 , a prominent ridged-ice 
formation on the subject aircraft wing was selected for 
detailed flow analysis using 2-dimensional(2-D), as well as, 
3-dimensional(3-D) Navier-Stokes computations. This con- 
figuration was selected because it caused the largest lift 
decrease and drag increase among all the ice shapes investi- 
gated in the earlier study. A grid sensitivity test was per- 
formed to find out the influence of grid spacing on the lift, 
drag, and associated angle-of-attack for the maximum lift 
(C, ). This study showed that grid resolution is important 

and a sensitivity analysis is an essential element ol the pro- 
cess in order to assure that the final solution is independent 
of the grid. 

The 2-D results suggested that a severe stability and control 
difficulty could have occurred at a slightly higher angle-of- 
attack(AOA) than the one recorded by the Flight Data 

Recorder (FDR) [2J . This stability and control problem was 
thought to have resulted from a decreased differential lift on 
the wings with respect to the normal loading for the configu- 
ration. The analysis also indicated that this stability and con- 
trol problem could have occurred whether or not natural ice 
shedding took place. 

Numerical results using an assumed 3-D ice shape showed 
an increase of the angle at which this phenomena occurred of 
about 4 degrees. As it occurred with the 2-D case, the trail- 
ing edge separation was observed but started only when the 
AOA was very close to the angle at which the maximum lilt 
occurred. 


Nomenclature 

MVD Median Volume Diameter (in p/n ) 

3 v 

LWC Liquid Water Content (in g/m ) 

CPS Control Point Smoothing 

CFD Computational Fluid Dynamics 

IGES Initial Graphics Exchange Specifications 
LW Left Wing 

RW Right Wing 

NTSB National Transportation Safety Board 

FDR Flight Data Recorder 

IRT Icing Research Tunnel (at NASA Lewis) 

AOA Angle-of-Attack, Degrees 

B-B Badwin-Barth Turbulence Model 

S-A Spalart-Allmaras Turbulence Model 

C j 2-Dimensional Lift Coefficient 

C, Maximum 2D Lift coefficient 

C L 3-Dimensional Lift Coefficient 

Cj 2-Dimensional Drag Coefficient 

C D 3-Dimensional Drag Coefficient 

y j Distance from the wing surface to the first grid 

point in normal direction (minimum wall spacing) 
v + Reynolds number based on the typical velocity 

and length scales for the turbulence 

Introduction 

Aircraft performance degradation due to ice contamination 
remains a concern within the aviation industry. Recent acci- 
dents and incidents have shown that undetected ice accretion 
or ineffective ice removal methods can lead to altered perfor- 
mance characteristics and sudden loss of stability and con- 
trol. with the potential for the most severe consequences* 31 . 
As part of a response to a request for technical assistance by 
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the NTSB, the NASA Lewis Icing Branch has used simula- 
tion methods to examine the possibility of ice contamination 
being a contributing factor in an accident involving a com- 
mercial turbo-prop powered aircraft. This examination 
involved tests in the Icing Research Tunnel (IRT) at NASA 
Lewis to obtain the ice shapes which would have accreted 
under the flight conditions of the accident aircraft ^ fol- 
lowed by a numerical analysis effort to determine the possi- 
ble performance degradation associated with those ice 
shapes, which is the subject of this paper. Because of higher 
than normal turbulence levels, partly due to the existence of 
the heat exchanger and spray bars necessary for icing cloud 
generation, the IRT is not considered the best facility for 
post-ice-accretion aero-measurement. High fidelity aeroper- 
formanee data is thus obtained from other wind tunnel facili- 
ties. However, due to constraints of cost and of available 
time associated with model fabrication and for extensive 
wind tunnel tests. Computational Fluid Dynamics (CFD) 
was considered an alternative means of providing post-ice- 
accretion aeroperformance analysis. 

An earlier numerical study by the authors 111 , aimed at build- 
ing a foundation for flow analysis of iced airfoils, attempted 
to define the iced surface as closely as possible to the origi- 
nal geometry in order to make the numerical grid generation 
and simulation process accurate and efficient. That study 
showed that, for similar ice accretion times, an ice formation 
with a prominent ridge caused the most severe performance 
degradation compared to other ice shapes having only sharp 
bumps of similar sizes. Based on this information and the 
flight conditions recorded by the FDR^ 2 l a numerical study 
was performed to determine if the degraded aerodynamics 
resulting from such an ice contamination could have led to 
the type of control upset attributed to the accident aircraft. 

Approach 

Lee Shape Mode line and Grid Generation 

The first step in the study was the modeling of the ice shapes 
obtained in the IRT^. An ice shape with a ridge on the 
upper surface near the leading edge and a number of small 
bumps on the lower surface (Fig. la & lb) was obtained in 
the IRT at the mid-section of a vertically mounted (Fig lc) 
wing. At the centerline of this 6 foot test model, the airfoil 
was roughly a NACA 23015 and the chord length was about 
68 inches. The icing spray conditions which produced this 
ice shape were a median volume diameter (MVD) of 20pm , 
a total temperature of 26 °F , an angle-of-attack of 5° , a liq- 
uid water content (LWC) of 0.8g/m^and a 5 minute spray 
time. 


Among the ice shapes numerically tested, the prominent 
leading edge ridged-ice caused the most significant perfor- 
mance degradation. The height of the ridge is less than 1% of 
the chc rd (0.0074c). An ice shape of this size would have 
been hard to see by the pilots in the cockpit. The geometry 
shown n fig. la has several sharp corners and high curvature 
segments. Curve discretization and eventual generation of a 
quality field grid on this geometry were time consuming. A 
mathematical model for systematic surface smoothing was 
presented in reference 1. This approach was implemented 
in an interactive code, TURBO-GRD [5] to generate surface 
shapes with different smoothing levels. It constructs a 
smooth curve whose shape is controlled by a piece-wise lin- 
ear curve formed by selected discretized points. These points 
are called control points (CPs) as they control the shape of 
the curve they construct. A brief summary of this process is 
presented here. 

1 ) The digitized ice shape data is first read in and one CP is 
assigned to each digitized point. 

2) A curve is constructed using these CPs. The discretized 
data points are moved onto the new smoothed curve deter- 
mined ly the CPs. The points are then redistributed at 
equally spaced intervals along the curve. The number of dis- 
cretized data points on the iced segment is unchanged by this 
process We will call this a baseline curve or a curve with 
100% Control Point Smoothing (CPS). 

3) Using this baseline curve as a starting point, the number 
(or percentage) of CPs can be reduced, thus generating 
curves with various levels of smoothness. During this pro- 
cess, the shape and the number of points in the un-iced areas 
do not change. 

The current study showed that a 50% level of CPS or higher 
is required to adequately represent the ice shapes as mea- 
sured by having marginal influence (i.e. less than 5% varia- 
tion) on the resulting lift and drag values. All grid generation 
was per ormed using the commercial code GRIDGEN^ 7 ^ 

The 2-1 > grids used in this study for modeling the region 
around the iced airfoil are composed of two blocks where the 
inner and outer blocks have an overlapping interface 
between ().5c and 0.6c (c=chord) from the airfoil surface. 
The inner block had a C-type grid with a wake cut down- 
stream c f the trailing edge and a much denser distribution of 
grid points than the outer block. Both blocks had down- 
stream toundary set at 15.0c from the leading edge. (Fig. 2a 
and 2b) The C-type grid was also used for the outer block 
with the farfield boundaries placed at a distance of 15 chord 
lengths rom the body surface in all directions. This method 
of constructing 2-block grids allowed for easier control of 
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the grid generation process and better quality ot the resulting 
grid, especially near the complex iced surfaces. 

Flow solver and boundary conditions 

A general purpose Reynolds averaged Navier-Stokes flow 
analysis code, NPARC^ was used for the simulation. In 
NPARC, the Navier-Stokes equations are formulated as cen- 
tral-difference approximations with added artificial dissipa- 
tion and are solved using an approximate-factorization 
scheme which results in a scalar penta-diagonal matrix for 
steady state computations. Complex geometries can gener- 
ally be handled with ease by the multi-block capability and 
modular boundary conditions. Inviscid, laminar, and turbu- 
lent flows can be simulated for 2-D (or axi-symmetric) and 
3-D geometries. A capability to calculate the lift and drag 
was added to a subroutine for this analysis. The code also 
has Runge-Kutta and implicit subiteration schemes tor time 
accurate computations. For the simulation ot turbulent flows, 
NPARC offers algebraic, one-equation, and two-equation 
turbulence models. In this study, both the Spalart-Allma- 
ras l8] and the Baldwin-Barth [9] one-equation turbulence 
models were used. At the tar field boundary, a non-reflect- 
ing type boundary condition was applied. 


Grid Sensitivity Test 

In the abscence of detailed measurement data ot the pressure 
or velocity fields, grid sensitivity tests can be used to deter- 
mine the optimum grid density. In order to develop the high- 
est quality simulation of aerodynamic properties of interest, 
such as lift. C, and drag, a series of grids (for the inner 

‘win 

block) having different resolutions in both normal and 
stream-wise (circumferential) directions were constructed 
(the outer block was fixed with the dimension of 1 15 x 20). 
The NPARC code was then run, using these grids, and the 
aerodynamic quantities of interest were compared. For this 
grid sensitivity test, the Spalart-Allmaras turbulence model 
was chosen for its known robustness in airfoil/wing calcula- 
tions. Examples of the importance of grid sensitivity testing 
and how grid properties can affect the simulation results, 
especially at high AOA. can be found in Ret. 10 and 1 1 . 

The first set of grids constructed (Table 1) was used for an 
investigation of the normal direction sensitivity, which was 
followed by a study on the effect of minimum wall spacing 
(y 1 ) (Table 2) and finally by a study of the effect of packing 
grid points in the stream-wise direction (Table 3). In tables 1 
and 3, ‘s' represents the stream-wise direction and V repre- 
sents the normal direction respectively. 


Figure 3a shows that the lift values obtained using the sin 1 
grid varied by over 10% from those obtained using the sln2 
grid. On the other hand, the lift values changed only 0.26 - 
2.21% when it was refined to the level of sln3 grid (detailed 
numerical values can be found in table 4). A similar trend 
was observed for drag values as indicated in Fig. 3b and 
table 4. This suggests that the nisi grid did not have a suffi- 
cient number of points to predict the maximum lift value 
while showing that any further refinement beyond the level 
of the sln2 grid was not necessary. This led to a further 
investigation of the effect of minimum wall spacing on aero- 
dynamic performance parameters (see Table 2). 

Another purpose of the second study was to find an appropri- 
ate value of yl to be used for 3-D grid generation which 
would allow efficient and fast convergence while not sacri- 
ficing accuracy. Throughout this study, the computations 
were performed until the L2-residual dropped at least 3 to 4 
orders of magnitude and the lift value changed by less than 

io- 7 . 

In this second step of the analysis, all aspects of the grid 
except the value of yl were fixed. Using the wl grid as a 
baseline, the yl value was either decreased or increased to 
investigate the effect on the lift. Figure 4 and Table 5 show 
that the increase of y 1 from 2.0 x JO' 6 to 5.0 x I O' 6 resulted in 
0.27 - 1.36% change in the lift value below 9 degree AOA 
but resulted in a 4.53% increase at 1 1 degree AOA. In the 
case of decreased yl (w3 grid), a simitar trend was noted 
except that it resulted in a reduction of the lift by 0.29 - 
2. 1 6 7c up to an AOA of 9 degrees and a drop of 6.9 1 % at 11 
degree AOA. An investigation of the v + values showed that 

the average value of v + was approximately 2.7 for the w2 
grid, 1.0 for the wl grid and less than 1.0 for the w3 grid 
respectively. This result shows that the AOA for the maxi- 
mum lift was predicted for all three grids as 9 degrees and 
that some difference existed in the computations at the 1 1 
degree AOA. From this study, it was decided to use a yl of 
2.0 x JO ’ 6 for the 2-D studies. For the 3-D study, a y 1 of 5.0 
x JO' 6 was used for the grid generation to minimize CPU 
time. 

To draw the final conclusion on the choice of proper grid res- 
olution, a stream-wise (circumferential) direction grid sensi- 
tivity test was performed using the three grids listed in table 
3. In this case, only the stream-wise point density was 
changed to determine its effect on the aerodynamic perfor- 
mance parameters. 

As Fig. 5 and table 6 show clearly, there was less than a 1% 
change in the lift values regardless of the resolution. 
Changes in the drag were less than 1% except at the higher 
AOAs for the s3n2 grid ( 1 ,58% maximum). From this test, 


3 

American Institute of Aeronautics and Astronautics 



AIAA 99-0375 


we decided to use the dimension of 350 x 50 with the yl 
value of 2.0 x I0' 6 as the base line. 

IzdimmsiomL ends used for the airfoil with rideed-ice and 
aileron deflection 

According to the FDR, at the time of the NASA-defined con- 
trol upset, the left aileron deflection was 2.56 degrees down 
and the right aileron deflection was 2.74 degrees up. This 
configuration was intended for a right banking movement. 
The approximate angle-of-attack and Reynolds number were 
7.8 degrees and 10 million (based on the chord) respectively. 
Another greater aileron deflection that occurred after the 
defined control upset was 7.94 degrees aileron down on the 
left wing (LW) and 8.26 degrees up on the right wing (RW). 

Four grids for the airfoil with ridged-ice accretion, reflecting 
these four different angles of aileron deflection, were gener- 
ated using the guidelines developed from the grid sensitivity 
tests. More grid points were used near the aileron in the 
stream-wise direction over the grids used in the sensitivity 
study. Another set of four grids for airfoils without ice but 
having the same aileron deflection angles were also gener- 
ated for comparison to the iced airfoils. 

3 -dimensi onal arid veneration 

An IGES-type surface data of the turbo-prop aircraft is pic- 
tured in a shaded mode in Fig 6a. The geometry data was 
obtained from the manufacturer and some assumptions were 
made to generate the grids for both iced and un-iced wings. 

A number of assumptions are worth noting: First, the avail- 
able data did not have information about the wing tip geome- 
try. Figures 7c-d show the approximated wing tip used in the 
analysis. Second, as shown in Fig 6b, only the outboard sec- 
tion of the wing (approximately 5.4 m in the span-wise 
direction) was used for the current analysis. Third, the 
downward deflected aileron on the left wing (only 2.56 
degree deflection angle) was modeled as having no gap with 
the wing surface in the span-wise direction. This was accom- 
plished by creating a smoothly connected surface shape 
shown in Fig 7e (see also Fig 7c). The idea behind this was 
that the existence of a fence on the real aircraft between the 
aileron and the outboard flap next to it would have prevented 
any span-wise flow. Fourth, for the modeling of an assumed 
3-D ice shape, the height and the shape of the ridged-ice on 
the upper surface and of the bumps on the lower surface were 
kept constant (Fig 7a & 7b). They were projected in both 
span-wise directions from the location where the ice shape 
data was taken in the IRT (see fig 6b). Since the wing did not 
have a large sweep angle and the maximum thickness was 
larger near the root, the root area had a lower collection effi- 


ciency. Thus, it was considered to be a plausible assumption 
to use i constant cross-section ice shape and to smoothly 
taper it near the wing tip (Fig. 7d). 

Figure 7f shows an example of grid lines around the ridge 
and a typical velocity profile showing separation behind it. 
Based on the above assumptions, 3-D grids consisting of 
approximately 400,000 and 1 million grid points in a single- 
block format were generated for the clean(un-iced) and iced 
wings with an aileron deflection angle of 2.56 degrees. The 
symmetry plane was located at z = 4.5 m (z=0 being the 
fuselage centerline) and the span-wise location of the aileron 
was between z = 6.3 m and 9.6 m. The 3-D computational 
results were compared to the 2-D results at the z=6.59 m 
location where the 2-D IRT ice shape was taken. 


Discussion of Results 


Aileron defl ection # ! (2.56 down LW \ 2.74 up RW) 

The figures Fig. 8a - 9b show the lift and drag coefficients 
versus angle-of-attack obtained by steady state computations 
using two different one-equation turbulence models, Spalart- 
Allmara; (S-A) t8] and Baldwin-Barth (B-B) 191 . This calcula- 
tion was performed to investigate differences in the numeri- 
cal prediction resulting from the application of these two 
models. Both figure 8a(S-A) and 9a(B-B) show that the ice 
caused a slight lift decrease at low AOAs and a further sig- 
nificant decrease at higher angles. 

The open circles represent lift values for the un-iced LW 
with 2.56 degree aileron down and the squares represent 
those ol un-iced RW with 2.74 degree aileron up. We 
assume ihat no ice shedding occurred on any of the wings. 
The lift was computed every 2 degrees between 3 and 15 
degrees for the un-iced wings. Without the ice, as the lift 
curves show, the airplane would not have had any problem 
banking :o the right with the higher lift on the LW. The iced 
wing computations were done for AOA of 3, 5, 7, 7.8, 9, 10, 
and 1 1 cegrees (7.8 degree was the AOA where the control 
problem occurred according to the FDR). 

A close examination of the lift change on the iced LW pre- 
dicted b> the S-A turbulence model in Fig 8a shows that the 
peak value occurred around 9 degrees. The maximum lift 
was predicted to occur at a slightly lower AOA, (7.8 degrees: 
Fig 9a), vhen the B-B model was used. The maximum lift 
on the iced RW occurred at approximately 10 degrees for 
the S-A node! (Fig 8a) and 9 degrees for the B-B model (Fig 
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9a). Also noticeable from the two graphs is that a reversal 
of lift values occur at around 1 1 degrees and 10.5 degrees 
respectively. The lift values on each iced-wing become very 
close at around 10 - 10.5 degree AOA and it could have 
caused difficulty in aileron control. 

Now. if we consider the possibility of ice completely shed- 
ding off the RW. as was indicated by the experiment in the 
IRT* 41 . the reversal of the lift values could occur at a lower 
angle-of-attack. In this case, 9 (S-A) or 8 (B-B) degrees 
would be the approximate angles where the reversal would 
have happened. These are angles fairly close to those 
recorded by the FDR. Computation for higher than 1 1 degree 
angle-of-attack was not attempted because of the expected 
unsteady nature of the flow. 

Figures 8b and 9b show the drag coefficient prediction using 
the two turbulence models. They indicate a considerable 
increase in the drag but show no cross-over of the drag val- 
ues for the iced wings. Both figures indicate that the LW 
showed consistently higher drag values than the RW. This 
supports the FDR record that there was a tendency for yaw in 
the counter-clockwise direction before the accident hap- 
pened. A little closer examination of the above four graphs 
also show that the S-A model predicted slightly higher lift 
and lower drag values than the B-B model. Though different 

f 1 21 

from the current study, a previous computational study 
using these turbulence models suggests that difference in the 
spreading rates (or amount of mixing) in the wake regions 
could be a contributing factor to these differences. Figures 
lOa-c show the leading edge and trailing edge separation 
predicted by the S-A turbulence model. The trailing edge 
separation starts at an AOA ol 7 degrees and this separation 
region grows as the AOA increases. When the AOA reaches 
9 degrees, the trailing edge separation covers almost 70% of 
the upper surface. Figure 1 1 shows the turbulent viscosity 
contour at the 9 degree AOA. 


Ailenm deflection # 2 [down UM. UAL up 8.26 RW1 

A higher angle aileron input was applied to attempt to bank 
to the right after the defined control upset. This condition 
was also simulated using the 2-D Navier-Stokes analysis 
procedure (Fig 12 - 15b). Unlike the above lower angle set- 
ting, the cross-over of lift was not observed, but, a consider- 
able decrease of the lift differential was predicted (Fig 14a 
and 15a). As it occurred with setting #1, shedding of the ice 
from the RW would have still caused the lift reversal, at a 
slightly higher AOA of 10 degrees. However, under aileron 
setting #2, the effectiveness of the aileron might have been 
questionable even at lower AOA again raising the possibility 
of the control dilemma. The S-A turbulence model predicted 


higher lift and lower drag values than the B-B model with 
the exception of the 1 1 degree AOA case. At this angle, the 
S-A model predicted a higher value for both lift and drag 
than the B-B model. 


7-D flow analysis 

Due to the time requirements for the calculations, the 3-D 
analysis was limited to only one downward aileron deflec- 
tion angle of 2.56 degrees on the left wing for both iced and 
clean cases. This calculation was performed to investigate 
whether there were any 3-D effects which altered the 2-D 
flow characteristics. The S-A model was used for all 3-D 
computations and the results were compared to those of the 
2-D computations using the same S-A model. 

The 3-D computation showed that the a f occurs at about 

1 3 degrees, approximately 4 degrees higher than the 2-D 
case. At 9 degree AOA, where the 2-D calculation showed 
considerable change in the pressure coefficient plot on the 
upper surface (Fig 16a), the 3-D calculation showed little 
difference between the clean and the iced wings (Fig 16b). 
The comparison of the 2-D and 3-D pressure coefficient 
shows some differences just aft of the ridge and at about 0.6 
x/c on the upper surface (Figure 16c). The rake profile and 
Mach contour plots from the 3-D analysis at an AOA of 1 3 
degrees (Fig 17a & 17b) indicate that no large scale separa- 
tion occurs near the trailing edge on the upper surface. At 1 5 
degree AOA though, a trailing edge separation starts, as 
shown in the next series of graphs depicting pressure coeffi- 
cient (Fig 18), span-wise variation of the pressure along 
lines parallel to the trailing edge (Fig 19a & 19b). Mach con- 
tour (Fig 20a) and velocity profile (Fig 20b). The rake profile 
at 15 degrees AOA (Fig 21) shows a large separation region 
on the upper surface as well as circular motion of the fluid 
parallel to the wing surface. This occurs near the span-wise 
region where the aileron was deflected. The change of span- 
wise pressure distribution shown in figure 22 indicates that a 
small scale trailing edge separation started at an AOA of 1 3 
degrees and intensified at 15 degree AOA. The differences in 
the lift and in the drag between the clean 3-D wing and the 
iced 3-D wing are shown in figures 23 and 24. The 3-D lift 
and drag are defined as 

^ L H r - D 

Cl =— andC„- — 

where q ^ is dynamic pressure and S represents a reference 

area (wing span x mean chord) and L and D are total lift 
and drag of the entire wing respectively. The causes of the 
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difference between the 2-D and 3-D computations are proba- 
bly due to one or more of the following. 

1 ) Relatively poor resolution of the 3-D grids. 

2) Proportionally smaller ridge height at the inboard section 
ol the wing could prevent the trailing edge separation at the 
lower AOA. Due to the 3-D nature of the How, this could 
affect the flow in the outboard section as well. 

3) The smoothly connected surface constructed for the 
deflected aileron could prevent possible 3-D effects in the 
span- wise direction. 

4) Downwash caused by the tip vortices increase the adverse 
pressure gradient near the leading edge (Fig 16c) producing 
an effect analogous to induced camber* 1 *^. 

Concluding Remarks 

The purpose of this study was to perform a post-ice-accre- 
tion CFD analysis of the contaminated airfoil/wing surfaces 
of a turbo-prop aircraft under the reconstructed icing condi- 
tions from an accident. This analysis was performed in order 
to obtain some qualitative trends and to provide insight into 
the aerodynamics that may have led to a control upset. The 
grid sensitivity tests which preceded the numerical simula- 
tion showed that the prediction of accurate lift and drag val- 
ues as well as AOA of maximum lift can be affected by the 
grid resolution. The 2-D analysis indicated that the control 
upset could have occurred with or without the complete ice 
shedding at or slightly higher than the angle-of-attack 
recorded by the FDR. The performance degradation was 
observed to be a result of a combination of trailing and lead- 
ing edge separation. In the case of 3-D analysis, the trailing 
edge separation was observed to start near the maximum lift 
point. The 3-D results also showed that the maximum lift 
occurred 4 degrees later than the 2-D case. More work is 
needed in 3-D ice shape modeling and in grid refinement to 
understand the differences between the 2-D and 3-D results. 
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Table 5. Grid sensitivity test by varying the minimum wall spacing 
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Table 6. Effect ol streamwise direction grid refinement 
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Fig. lb Rough ice shape on the lower side of the Fig. lc Vertically mounted wing in the test section 

turboprop wing - IRT experiment 



Fig. 2a Two-block grid system used for 2-D analysis Fig. 2b Detailed grid near the ridged-ice 
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Fig. 3a Grid sensitivity test in normal direction 

Ridged-ice, S-A turb. model, inner block grid density changes 


Fi 3b Grid Sensitivity in Normal Direction 

Ridged-ice, S-A turb. model, inner block grid density changes 




G — G slnl grid- 350 x 35 
□ — 0 si n2 grid * 350 x 50 
A — A si n3 grid - 350 x 70 


0 — G si nl grid - 350 x 35 
Q — 0 Sln2 grid - 350x50 
A — A si n3 grid - 350 x 70 
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: ig. 4 Grid sensitivity test using minimum wall spacing 

Various y 1 , S-A turb. model used 

i.io i — ~ ^ i > ’ » ■ — ■ 


Fig. 5 Streamwise direction grid sensitivity test 

Ridged-ice, S-A turb. model, inner block grid density changes 


G — G wl grid, yl = 2e-6 
Q — 0 w2 grid, yl = 5e-6 
A — A w3 grid, yl = 1e-6 


G — G s1n2 grid - 350 x 50 
B — 0 s2n2 grid - 391 x 50 
A — A s3n2 grid - 435 x 50 
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Fi«. 6a Overall picture of the turboprop aircraft used for the numerical analysis in a shaded mode 



Fig. 6b Three dimensional wing modeling showing the area for numerical analysis 
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Fig. 7a Ice shape at the LE of the wing Fig. 7b Ice shape modelling on the lower surface 



Fig. 7e Modelling ot the aileron at 2.56 deg. down Fig. 7f Grid system at the rideged-ice and the 

separation behind it 
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g. 9a Lift Comparison for airfoil with ridged ice, B-B turb. model 



Fig, 9b Drag Comparison for a 
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ig. 15a Lift Comparison for airfoil with ridged ice, B-B 

Aileron up 8.26 deg. on the RW. down 7.94 deg on 
I ^ ' ’ 
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Fig. 18 Cp change due to various AOA for 3-D iced wing 


Aileron down 2.56 deg LW, S-A turb model, at z = 6.59 m 
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